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Abstract 

A detailed study of orbital and radial excited states in D, Dg, B and 
Bg systems is performed. The chiral quark model provides the frame- 
work for the calculation of pseudoscalar meson (vr, A', ...) hadronic 
transitions among heavy- light excited and ground states. To calcu- 
late the excited states masses and wavefunctions, we must resort to 
a relativistic quark model. Our model includes the leading order cor- 
rections in l/m(^(.,b) (e.g. mixing). Numerical results for masses and 
light hadronic transition rates are compared to existing experimental 
data. The effective coupling of the chiral quark model can be deter- 
mined by comparing with independent results from lattice simulations 
{g\ = 0.53 ± 0.11) or fitting to known widths {g\ = 0.82 ± 0.09). 



1 



1 Introduction 



Although heavy quark spectroscopy is now a rather mature subject, a num- 
ber of interesting issues remain. In particular, the detailed properties of 
the excitation spectrum of heavy-light mesons {D, Ds, B, Bs) and their light 
hadronic transitions are yet to be fully understood. Experimentally, much 
of this excitation spectrum remains to be observed. Only the ground state 
S-waves and a few of the ji = 3/2 P-waves are presently well established. 
However, many of these states will be accessible in the present and future B 
factories: CLEO, BaBar, Belle, CDF, DO, BTeV and LHC B. Besides further- 
ing our understanding of QCD dynamics, the detailed study of these excited 
states may have practical benefits. For example, triggering on excited states 
may provide an efficient method of same side B tagging in hadron colliders. 
Tagging is essential to the study of CP violation in the B system. 

The theoretical tools available to determine the properties of excited 
states in heavy-light systems include Heavy Quark Effective Theory (HQET) 
and low energy chiral effective theory Unfortunately, these tools are not 
sufficient to determine the detailed properties of these states. Lattice gauge 
theory is the only existing technique which allows the systematic study of 
all the aspects of QCD in heavy-light systems. Detailed studies of the P- 
wave excited heavy-light states within the quenched approximation already 
exist 10]. Future lattice studies will provide more insight into the nature of 
QCD dynamics as well as the masses and static properties of heavy-light 
hadrons. 

It is clear that a model to estimate the hadronic transitions from excited 
state to ground states would also be very useful. Such a formalism has 
been developed and applied extensively to transitions in heavy-heavy {QQ') 
systems . However, heavy-light mesons are more difficult because the light 
quark is subject to the full nonperturbative QCD dynamics. One possible 
approach to providing a framework for these hadronic transitions is to use 
the chiral quark model [^]. This has been suggested by Goity and Roberts 

In this paper we closely follow the work done in Refs. [§, 0]- To compute 
the masses and wavefunctions of the excited states we use a Dirac equation 
for the light quark in the potential generated by the heavy quark (including 
first order corrections in the heavy quark expansion and mixing effects). 
We then use these masses and wavefunctions to compute the hadronic decay 
amplitudes of excited heavy mesons in the context of a chiral quark model Q . 

The main differences between the present and preceding works are in the 
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choice for the parameters of the chromoelectric potential and in the inclusion 
of mixing effects both in the spectrum and in the decay amplitudes. Moreover 
we use our results for the radial wavef unctions of the excited mesons to make 
a comparison with recent lattice results. From the comparison we extract 
an estimation for g\, the effective coupling of the quark to the pseudoscalar 
mesons. We find g\ = 0.53 ± 0.11 

We present numerical results for the low-lying spectrum (excited states 
up to the 3S states). We also compute the pseudoscalar meson hadronic 
transitions for these states as a function of the chiral quark model effective 
coupling constant. Comparing our results with recent experimental width 
measurements we estimate this effective coupling g\ = 0.82 ± 0.09. 

In Section 2, we discuss our determination of the spectrum of excited 
states. Our notation, the choice of the potential, inclusion of mixing and 
other order l/ruh corrections are explained. Details of the masses and wave- 
functions are presented for the low-lying excitation spectrum. Comparison is 
made with present experimental data. Our treatment of hadronic decays is 



described in Section 3. The analytic results are summarized in Eqs. (|3^|3^ 
Explicit expressions for the coupling coefficients appearing in these equations 
are given in Appendix A. Also in Section 3, details of the partial rates for the 
IS and IP states in the D, Ds, B and Bg systems are presented. Again com- 
parison is made with present experimental data. A complete list of remaining 
results for masses and partial decay widths is reported in Appendix B. 



2 Spectrum 

2.1 Basic Model and Notation 

The general Hamiltonian of the heavy-light system can be expanded in powers 
of {l/mh) 

n = 7^(0) + ^^(1) + 1 ^(2) + 

rrih m% 

However, even within the heavy quark limit, the general form of the zeroth 
order Hamiltonian, 7i°, still involves the full nonperturbative QCD dynamics 
for the remaining degrees of freedom (including light quark pair creation and 
gluonic degrees of freedom). At present it can not be solved analytically. We 
are forced to resort to use a relativistic potential model for 7i°. 
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We model the most general heavy-light meson (in the D, Dg, B, Bg fam- 
ily), H, as a bound state of a light quark (g) and a heavy quark (h). The 
heavy quark is treated as a static source of chromoelectric field and the only 
quantum number associated with it is its spin. The light quark is treated 
relativistically and its state is described by the wavefunction ipn,i,j,m{r, 9, ip). 
In analogy with the hydrogen atom, we introduce the following quantum 
numbers: 

• n, the number associated with the radial excitations; 

• the orbital angular momentum; 

• j, the total angular momentum of the light quark; 

• m, the component of j along the 'z axis; 

• J, the total angular momentum of the system; 

• M, the component of J along the 'z axis; 

• S, the spin of the heavy quark along the 'z axis; 

The parameters of our model are the masses of the light quarks (m^ for 
q = u,d or s), the masses of the heavy quarks {rrih for h = c or b) and the 
chromoelectric potential of the heavy quark {V{r)). 

The total wavefunction of the system can be decomposed as follows 

^n,£,j,.J,Mi^^ ^^V) = ^ C*/'^! 5^n,£j,m(r, 6*, if) ® (2) 

where C'^^'^i „ are the usual Clebsh-Gordan coefficients and fs is a two com- 

ponent spinor representing the heavy quark. Eq. is a solution of the 
following eigenvalue problem 

(3) 

where Ti is the Hamiltonian of the system. The energy levels in Eq. (|^) do 
not depend on M because of rotational invariance. 
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We rewrite Eq. (0) introducing the most general parameterization for the 
four spin components of the light quark wavefunction 



Aj,J,M 



C 



J,M 



m m- 



\ fn/,ji.^)^2j-l,j,m-m+i 



■ Y 

2j-e 



(4) 



Here Y^{9,Lp) are spherical harmonics that encode the angular dependence 
while fn/,jij)i fn,i,j{'^) real functions that encode that radial dependence. 
k'^- ^ and kj- ^ are fixed, up to an overall phase, by imposing a normalization 
condition. Our choice of the phase is such that 

2.2 Choice of the potential 

Within our basic framework, Ti'^'') is given by the relativistic Dirac Hamilto- 
nian 

7^(0) =70(-^^ + mg) +V(r) (6) 

and the rotational-invariant potential is the sum of a constant factor (Mh), 
a scalar part (K) and (the zeroth component of) a vector part (K) 

Vir) =Mh + 7°V;(r) + K(r) (7) 

The constant Mh is a an overall energy shift that depends on the heavy quark 
flavor and, in general, it is not equal to m^, as often assumed in the literature. 
For this reason we consider m/j and two independent parameters of the 
model. 

Asymptotic freedom suggests that at short distances the potential is dom- 
inated by a vector part that asymptotically approaches a Coulomb potential, 
V ~ ~ 1/r. On the other hand lattice simulations indicate that at 



k^ 
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large distances the potential is confining, scalar and asymptotically linear, 
V Vs ^ r. 

The naive assumption about a short distance Coulomb-like divergent be- 
havior of the potential is doomed because it gives rise to ultraviolet diver- 
gences (as discussed in Ref. P] and Ref. [§]). It this context the divergence 
arises in the 1 / rrih correction to the energy and it is due to the inconsistency 
of a static point-like source (the heavy quark) within a relativistic framework. 
One solution is assuming that the heavy quark is static but not point-like, 
therefore the potential that it generates is a convolution of the Coulomb-like 
potential and the square of the heavy quark wavefunction (peaked around 
the center of mass of the system and smeared within some small length scale 

More generally, one is allowed to cure this divergence by regulating the 
potential close to the origin (on a length scale of the order A~^). Different 
choices for the regulator are allowed and they do not affect the physics we 
want to describe, providing that is small enough. The values of the 
parameters that appear in the Hamiltonian, on the contrary, depend on this 
choice since they run with A. In fact, to obtain the same spectrum, different 
choices for the regulator imply different fitting parameters. 

We chose to regulate the vector potential by assuming a Gaussian shape 
for the wavefunction of the heavy quark, <l>(x) = exp(— x^A^/2), and with 
this choice 

Vv{r) = -\ I |$(x)|' y^d?x = -l%rf(Ar) (8) 
6 J \y — 6 r 

For the scalar potential we assume a simple linear form 

Vs{r)=br + c (9) 

We observe that c is not a physical parameter since it can be absorbed into 
the definition of m^. For this reason c will be omitted from now on. 
Summarizing, the nine parameters of our model are 

as,X,b,mu,ms,mc,Mc,mb,Mb (10) 

where = rrid and are mass parameters for the light u, d and s quarks 
respectively, equivalent to constituent quark masses shifted by the constant 
amount c of eq. (P), which is undetermined in our model, is the mass of 
the c quark with the corresponding energy shift. Analogously for the h 
quark. 
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2.3 1/m/i correction 

For any given set of input parameters we solve the eigenvalue problem, 
Eq. (H), using the Hamiltonian of Eq. @ and the potential specified by 
Eqs. (0), (H) and (|^). In this way we determine the radial wavefunctions 
and f^^^ j associated to the energy levels E^lj- We then compute l/ruh 
corrections to the energy levels in first order perturbation theory 



n,£,j,J 



E. 



(0) 



(11) 



where, ignoring for the moment the mixing of the states. 



SE. 



(1) 

n,e,j,J 



M 



X 



'12) 



The analytical expression for H.^^^ has been derived in Ref. using the 
Bethe-Salpeter formalism. In terms of the radial wavefunctions (after the 
analytical integration of the angular part), we rewrite SE^^^ as a sum of 
three contributions 



J, J 



(13) 



These terms are respectively: 
• The kinetic energy 



A 



n,£,j,J 







fn,e,j 



dl + -dr 



dl + -dr 



+ 



+ 



f 
J 1 







(14) 



fl,^^ rMr (15) 



with £ = 2j - 

A shift due to spin-orbit interaction 



B. 



n,e,j,J 



fn,e,j ( - ) fn,e,j fn,e,j \ 



+ 2 



fnl.j 



r dr 



(16) 
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for j = £ + |, or 

Bn,e,j,J fn,i,j (^^'■ + ~^^ fn,i,j-fn,i,j (^^r 



(17) 



for J =£ - |. 



• The hyperfine splitting 

Cn,i,,J = (-1)'"' j^m.)fl,jU, rMr (18) 

2.4 Mixing 

In Eqs. ([TT|) and (|12D we assumed that the Hamiltonian was diagonal. This is 
not the case because the — interaction term in the Hamiltonian mixes states. 
In general correction terms can mix any states with the same total angular 
momentum, J, and parity, P, However, there are only two types of sizable 
mixings. Large mixing can occur for pairs of states ^n,i,j,j,M and '^n',i',j',j,M 
with: (1) n = n', i = i' and J + 1 = j' ~ | = i = J (i.e. mixing within a given 
n and / multiplet); or (2) n + 1 = n', £ + 2 = f and j + ^=j'-^=i = J 
(e.g. S-D mixing). The off-diagonal term in the Hamiltonian that mixes such 
pairs of states has the form 

/ *l,..,J,M(^)^^'^^n',r,,',J,M(x)d^X (19) 



1 VJ(J+1) 



' ' nih 2J+1 

and it induces a mixing in the wavefunctions and in the energy levels 



'4^n,e,j,m 




phys. ^ 








phys. ^ 


En',e',j',J y 






En',i',j',j)/'^ 



"2K -'- 7 \ '^n',e',j',m 



+ 0{e^) (20) 



^n,^j,J + 2A^^ \ ^ ^^^2) (21) 
En',e',j',J - 2A / 



We find that the effect of mixing is generally negligible except for the P 
waves, where the mixing among wavefunctions can be of the order of 10%. 
In Tables |l]-| we report the value of = ^e/A for each excited state. It 
measures, in percent, the contribution of the mixing to the wavefunction. 
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TT 

H 


wLj) 


^exp. 


ti 




(/)(%) 


D 


(la^o) 


1.865 


1.895 


1.868 




D 


/I i o \ 


2.007 


1.895 


2.005 




D 






2.282 


2.377 




D 


(I2P1) 


2.422 


2.253 


2.417 


-10.92 


D 




2.459 


2.253 


2.460 




D 


(I2P1) 




2.282 


2.490 


10.92 


D 


(22 5*0) 




2.447 


2.589 




D 


(22 5l) 




2.447 


2.692 


2.17 


D 


(12-D2) 




2.504 


2.775 


-5.41 


D 


(12-D1) 




2.553 


2.795 


-2.17 


L> 






n /I 
Z.0U4 


z. lyy 




D 


(11^2) 




2.553 


2.833 


5.41 


D 


(25Po) 




2.683 


2.949 




D 


(2iPi) 




2.679 


2.995 


-10.70 


D 


(2iP2) 




2.679 


3.035 


1.79 


D 


(2^ Pi 




2.683 


3.045 


10.70 


D 


(l^i^s) 




2.709 


3.074 


-3.17 


D 


(I5P4) 




2.709 


3.091 




D 


(I5P2) 




2.760 


3.101 


-1.79 


D 


(iji^a) 




2.760 


3.123 


3.17 


D 


(3 2 5*0) 




2.823 


3.141 




D 


(32 S*!) 




2.823 


3.226 





Tabic 1: Tabulated spectrum for D mesons. E'^ denotes the lowest order 
energies. E^^^^' includes all the order l/rrih corrections. (All units in GeV). 
The mixings between lowest order states are denoted by 0. 



9 



H ( 


n^Lj) 


^cxp. 


^0 


^phys. 




Ds 


(la^o) 


1.969 


1.988 


1.965 




Ds 


(l^^l) 


2.112 


1.988 


2.113 




Ds 


(l^Po) 




2.374 


2.487 




Ds 


(ifPl) 


2.535 


2.353 


2.535 


-11.62 


Ds 


(itPa) 


2.573 


2.353 


2.581 




Ds 


(l^pl) 




2.374 


2.605 


11.62 


Ds 


(2^5o) 




2.540 


2.700 




Ds 


(2^Si) 




2.540 


2.806 


1.97 


Ds 






2.606 


2.900 


-6.11 


Ds 


(liDi) 




2.648 


2.913 


-1.97 


Ds 


H^D,) 




2.606 


2.925 




Ds 


(liDa) 




2.648 


2.953 


6.11 


Ds 


(2^Po) 




2.777 


3.067 




Ds 


(2iPi) 




2.775 


3.114 


-10.58 


Ds 


(2tP2) 




2.775 


3.157 


1.81 


Ds 


(2^Pi) 




2.777 


3.165 


10.58 


Ds 






2.812 


3.203 


-3.60 


Ds 


(liPi) 




2.812 


3.220 




Ds 


(1§P2) 




2.857 


3.224 


-1.81 


Ds 


(iji^s) 




2.857 


3.247 


3.60 


Ds 


(32 So) 




2.917 


3.259 




Ds 


(32 5*1) 




2.917 


3.345 





Table 2: Tabulated spectrum for Ds mesons. (All units in GeV). Notation 
as in Table |l|. 
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H 












B 




5.279 


5.288 


5.279 




B 




5.325 


5.288 


5.324 




B 


(liPi) 




5.646 


5.700 


-6.00 


B 


(l^Po) 




5.675 


5.706 




B 


(1^P2) 




5.646 


5.714 




B 


(l5Pi) 




5.675 


5.742 


6.00 


B 


(2^5o) 




5.840 


5.886 




B 


(255i) 




5.840 


5.920 


0.69 


B 


{1ID2) 




5.897 


5.985 


-1.96 


B 






5.897 


5.993 




B 


(iIdi) 




5.946 


6.025 


-0.69 


B 


{1ID2) 




5.946 


6.037 


1.96 


B 


(25Po) 




6.076 


6.163 




B 


(2iPi) 




6.072 


6.175 


-9.11 


B 


(2iP2) 




6.072 


6.188 


0.50 


B 


(25Pi) 




6.076 


6.194 


9.11 


B 






6.102 


6.220 


-0.99 


B 


(liF4) 




6.102 


6.226 




B 


(1^P2) 




6.153 


6.264 


-0.50 


B 


(iji^s) 




6.153 


6.271 


0.99 


B 


(SaS'o) 




6.216 


6.320 




B 


{32 Si) 




6.216 


6.347 





Table 3: Tabulated spectrum for B mesons. (All units in GeV). Notation as 
in Table |l[ 
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H 














Bs 


(1^5o) 


5.369 


5 


381 


5.373 




Bs 


(i^^i) 


5.417 


5 


381 


5.421 




Bs 


(l5Po) 




5 


767 


5.804 




Bs 


(liPl) 




5 


746 


5.805 


-7.19 


Bs 


(liP2) 




5 


746 


5.820 




Bs 


(l5Pi) 




5 


767 


5.842 


7.19 


Bs 


(25^o) 




5 


933 


5.985 




Bs 


(2i^i) 




5 


933 


6.019 


0.64 


Bs 


(11^2) 




5 


999 


6.095 


-2.31 


Bs 


(1ID3) 




5 


999 


6.103 




Bs 






6 


041 


6.127 


-0.64 


Bs 


(11^2) 




6 


041 


6.140 


2.31 


Bs 


(25Po) 




6 


170 


6.264 




Bs 


(2iPi) 




6 


168 


6.278 


-9.81 


Bs 


(2iP2) 




6 


168 


6.292 


0.51 


Bs 


(2lPi) 




6 


170 


6.296 


9.81 


Bs 


(I5P3) 




6 


205 


6.332 


-1.15 


Bs 


(I5P4) 




6 


205 


6.337 




Bs 


(ltP2) 




6 


250 


6.369 


-0.51 


Bs 


(1JP3) 




6 


250 


6.376 


1.15 


Bs 


(325'o) 




6 


310 


6.421 




Bs 


(3^5i) 




6 


310 


6.449 





Table 4: Tabulated spectrum for Bg mesons. (All units in GeV). Notation 
as in Table |l|. 
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B mesons B mesons 

s 

6.6 I , , , , 1 6.6 I , , r 




5.2 - - 5.2 - 



5.0 I ' ' ' ' 1 5.0 I ' ' ' L 

1 2 3 1 2 3 

L L 



D mesons D. mesons 




1 2 3 1 2 3 

L L 

Figure 1: Computed spectrum of excited states in the {D, Dg, B, Bg} family. 
The plot shows the spectrum before and after l/nih corrections (including 
mixing). These corrections are responsible for the hyperfine splitting. The 
horizontal axis is the orbital angular momentum of the meson {£). For each 
value of i and j there is a doublet of states {J — j — ^ and J — j + ^, with 
lower and higher energy respectively). 
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r(tm) 



r(fm) 



Figure 2: Radial wavefunctions for some excited states (for non-strange 
mesons). The continuum (dashed) hne refers to the /°(r) {f^{r)) function. 
These plots do not include the mixing contribution. 
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2.5 Determination of the parameters and predictions 



The nine parameters of our model, Eg. (p!OD , are determined numerically as 
follows: We define a function T of the input parameters that finds eigenvalues 
and eigenfunctions of Eq. (|]) using a fourth order Runge-Kutta formula, 
corrects the energy levels by including the 1/m/i perturbative corrections 
(including mixing effects) and returns 



X 



E 



observed states 



(22) 



where E^lj'j are the computed energy levels and mn,£j,j ± Sm are the mea- 
sured masses (with their experimental error) of the corresponding particles. 

We then minimize in its nine dimensional domain. We repeat this 
procedure with different sets of starting parameters until we are confident 
that we have found the absolute minimum. The experimental data used for 
the "observed states" in the fit are reported in the third column of Tables P 
Our best fit gives the following values for the parameters 



as 


0.339 


A 


2.823 GeV 


b 


0.257 GeV^ 


niu 


0.071 GeV 




0.216 GeV 


rric 


1.511 GeV 




1.292 GeV 


rrib 


4.655 GeV 


Mb 


4.685 GeV 



(23) 



The corresponding predicted spectrum is reported in the fifth column of 
Tables The best fit parameters reported here differ slightly from those 
reported in Ref. [1^ because we use here most recent value for mn/,j,j- 



We remark that m„ and are mass parameters which differ from the 
constuent quark masses for an overall undetermined constant shift. 

As a consistency check of our results we observe that the mass splitting 
rris — rriu — 140MeV comes out in agreement with naive expectations based 
on Gell-Mann-Okubo type relations. This difference also agrees, within 1%, 
with the corresponding splitting determined in ref. and used in ref. as 
input for their calculations. 
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-0.5 

X (fm) 




Figure 3: Orbitals for some excited B mesons. 



Remarkably mj, ~ with much better agreement than expected. More- 
over ~ 0.06fm is smaller than any other length scale involved in the 
problem, as was required. 

Figure ^ shows some of the computed radial wavefunctions for non-strange 
mesons, f^^i jir) and fn/ji''^)- These wavefunctions do not include l/rrih cor- 
rections, therefore they are the same for D and B mesons. The corresponding 
wavefunctions for strange mesons are very similar. 

Figure ^ shows density plots for each couple of independent spin compo- 
nents of some of the computed light-quark wavefunctions, |/°.^j^o^| in black 

2 

(blue) and fn,e,jYo''~ S^^y (red). They represent the analogous, in the 
heavy meson systems, of the orbitals of the hydrogen atom. 



2.6 Comparison with experiment 

The comparison of our results to the present experimental information on the 
excitation spectrum of the {D, Dg, B, Bg) mesons is given in Table ^ States 
which were used in the determining our best fit parameters are so indicated. 
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charmed meson masses (MeV) 



bottom meson masses (MeV) 



model exp 



Ref. 



model exp 



Ref. 



D 137 [a] 141,142 P 



B* - B 45 [a] 46 



D*-D 
Dl - D 
Di-D 



D* 



512 

622 596(53) |14| 
549 [a] 558(2) 



11 



D 592 [a] 594(2) PI 



Bo-B 
Bl - B 
Bi- B 
Bl- B 
B*-B 

{B** - B) 



427 
463 
421 
421 
435 



391(16) 

459(9) 

431(20) 

489(8) 

460(13) 

418(9) 



mm 

0[c] 



D' - 
D*' 



D 
-D 



721 
824 



772(6) PI 
not seen 



B' - 
B*' 



B 
- B 



607 
641 



16, 17 



Ds-D 97 [a] 99,104 
D*^ - D, 148 [a] 144 



m 
11 



Bs-B 
B* - B. 



94 [a] 90(2) |TT 
48 [a] 46 PI 



D 



sO 



512 



B 



Dl^-Ds 640 
Dsi - 570 [a] 



566(11 



Dl^ - Ds 616 [a] 605(2) 



-B, 

-B, 
-B, 
-Bs 

{Br - Bs 



sQ 

Bti 
Bsi 

^s2 



431 
469 
432 
447 



484(15) pl[c 



D'^ - Ds 735 
D*' - D, 841 



B's-Bs 
B*J - B, 



612 
646 



Table 5: The heavy-light spectrum compared to experiment. We report the 
difference between the excited state masses and the ground state [D or B) 
in each case. 

[a] Experimental input to model parameters fit. 

[b] Theoretical estimates for some of the mass splittings have been used as 
input. 

[c] Experimental signal is a sum over resonances with J = 0, 1, 2. 
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Our model is in excellent agreement with the better established P waves 
in the D and Dg systems. In particular the D\ (I2P1) fits the recent measure- 
ment of CLEO |T^. For the P- waves of the B meson systems the agreement 
with preliminary measurements is somewhat less impressive. 

However many of the existing experimental fits for individual masses 
of these states relied on patterns of masses for the ji = 1/2 states not 
found in our model. For example, our relativistic quark model predict 

1 i 1 > 1 3 I for a relatively broad range of parameters consistent with 

light spectroscopy^. Also, we obtain a splitting for the ^Pj states more than 
twice as big as the splitting for the aPj states. 

Preliminary results from L3 for the masses of P-wave B meson exci- 
tations are: 



mi 1 1 1 = (5.670 ± O.OlOstat ± 0.013syst)GeV (L3) (24) 
^2 ■ ^i,i,|,2 = (5-768 ± O.OOSstat ± 0.006syst)GeV (L3) (25) 

(26) 



m 



1,1,^,0 



The L3 results were derived using the constraint that m^-^^ 1 ^ 

]^ 3 2 — I 3 ^ = 12 MeV. This assumption is not realize(i in our model. 
Similar assumptions are needed for the extractions the masses of P-wave B 
meson excitations from OPAL |jl3| and CDF [Q. It would be interesting to 
reanalyze results using the pattern expected in this relativistic quark model. 



Finally, the observation of the D*' by DELPHI |15] is not consistent with 



searches by CLEO |16l and OPAL g 



2.7 Regge trajectories 

We find that En/j,j — Mh all lie on Regge trajectories parameterized by 
mJ + qni,j with m ~ 0.7 (as shown in Figure ^). This is a well understood 
phenomenon for light spectroscopy but, for light mesons, m ~ 27r6 (where h 
is the string tension). This is about a factor two bigger than we find. Our 
result can be explained in the non relativistic limit or, alternatively, in a 
simple and naive string picture: Ignoring the short distance behavior of the 
potential, we model the meson as a classical light quark attached through 
a string to the center of mass of the system. The energy of the system, E, 

"'^This result is known in the literature as spin-orbit inversion. It was first predicted by 
Schnitzer ||2^ and later by the models of Isgur and Ebert ||2^ . 
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Figure 4: Regge trajectories for some of our orbitally excited heavy mesons. 
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(i.e. the energy of the string) is related to the classical angular momentum, 
J, by -E^ = TcbJ. The factor two in the light meson picture can be explained 
with the fact that the latter rotate around a center of mass that is located 
at the middle of the string while for heavy-light mesons the center of mass 
coincides with one of the two ends of the string. 

This simple picture shows that the bulk of an heavy meson mass is domi- 
nated by the mass of the heavy quark plus the potential energy associated to 
the large distance interaction (oc br), and again gives support to our assump- 
tion that short distance behavior of the potential has a small contribution to 
the spectrum. 

3 Hadronic Transitions 
3.1 Transition amplitudes 

We start by considering the most general hadronic transition of the form 

H' ^H + x (27) 

where H' and H are two heavy-light mesons containing the same heavy quark, 
with wavefunctions '^n',e',j',j',M' and '^n/,j,j,M respectively and x can be any 
light meson with momentum p. Although we keep our formalism general, in 
this paper we only compute numerically decays in which x is a pseudoscalar 
meson belonging to the flavor octet (vr, rf). 

In the context of the chiral quark model this transition is mediated by 
an effective interaction of the form 

C,^, = -A_^^XM'^q, + Oid') (28) 

where fx can be identified with ~ 130MeV and g\ is an effective coupling. 
X = is the spin structure associated to the transition, i,j are 5'[/(3)flavor 
indices and 



M = V2 



(29) 



is the usual SU{?>)l^r invariant representation of the pseudoscalar mesons. 
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H' - 


^ H + x 


c 


iJ' - 


H + x 


c 


50 - 




1 


Bs- 


^B, + 7]s 


-2/VS 






1 


Bs- 


^ B^ + K 


V2 




-> SO + 7r± 


V2 


Bs- 


^ B^ + K^ 


V2 


50 - 




V2 


B'- 


-^Bs + K 


V2 






1/v^ 


- 


-^Bs + K^ 


V2 






1/v^ 


B,- 


^ B,, + r]i 


\/f + ^(i) 



Table 6: List of decay channels for B (or D) mesons with the corresponding 
flavor factor (. 



The transition mediated by this Lagrangian is associated to the following 
matrix element 

/^'^^(p) = ^ J j,M(^)Xe^P-vl>„, ,, ., ^, ^,,(^)d3;^ (30) 

where C is a coefficient that characterize the flavor structure of the decay. A 
list of all possible cases has been derived from Eq. (^) and is reported in 
Table The physical 77 and rj' are of course mixtures of the ideal ?78 and rji. 
In particular, 77 = ?78Cos(6'p) — r/i sin(6'p) where 9p ^ — lO.ldeg with a large 
uncertainty. In this work we ignore this mixing and we assume 6p = 0. The 
corrective multiplicative factor for a different choice for 6p can be derived by 
the reader using Table |. 

The exponential in Eq. (|30[) can be expanded in products of spherical 
harmonics and spherical Bessel functions, thus giving 

/-'-(p) = y: iir(p)<^:,^^<-(x,p) (31) 

Eq. ( prp implicitly defines the transition amplitude, Af^^^^{X,p), for a x 
in a given eigenstate of its angular momentum. By projecting the matrix 
X on the basis presented in the Appendix A, the transition amplitude can 
be rewritten as a linear combination of terms, each factorized into a radial 
part and a spin dependent part 



(32) 
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The coefficients c"^'^(X) depend on the quantum numbers of the mother 
and the daughter heavy mesons. Their exphcit expression is given in the 
Appendix A. The integrals are computed numerically. 

One can extend our analysis for octet pseudoscalar transitions to the 
approximately flavor singlet v( . In the large Nc limit the r/i combines with the 
octet to form a nonet. In that case the effective interaction in the Lagrangian 
takes the form 

A„t = ^g^Xr/ig, + 0{d^) (33) 

with X = and ~ (7^. This symmetry is badly broken in QCD. However 
it is reasonable to assume that in these transitions that the form (Eq. (p3|)) 
still holds. If one further assumes that the spatial wavefunctions of vr and 77' 
are the approximately equal, one obtains jr^' — fn- Hence the coefficient ( 
can be set to a/2/3 both for heavy strange and non strange decaying heavy 
mesons. 

The situation for decays in which x is a light vector mesons (p, uj, K*) is 
different. When compared to the pseudoscalar mesons, they have a different 
spin coupling to the quarks {X = / where is the polarization vector of the 
meson), a different effective coupling {gv 7^ Qa), and a different wavefunction 
ifp f-rr)- With these replacements Eq. ( ]32D remains valid for decays with 
emission of light vector mesons. [| The detailed study of these vector meson 
transitions is deferred to a future paper. 



3.2 Partial widths 



The partial width for the transition in Eq. (^) (for a light meson x emitted 
with total momentum p and angular momentum i^) is given by 



p 2 J + 1 rriH 

87r2 2 J' + 1 TTlH' 



(34) 



where m h' and ttih are the masses of the mother and daughter heavy mesons 
respectively. 



^It is possible to relate these the pseudoscalar and vector couplings and coefficients 
within the context of an approximate SU(6)w symmetry for the low lying states. 
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channel 


p. (MeV) rj((7i)2(keV) 


^D+(l^S'o) + 7rO 


1 42.8 ±.2 62 ±1 
1 39.4 ±.5 97 ±3 
1 38.1 ±.3 44 ±1 



Table 7: List of partial decay rates for IS D mesons. The measured D* and 



D masses [O] are used in these results. 



The total hadronic decay width (via a pseudoscalar meson transition) is 
defined simply as the sum of the partial widths: 



■^H' 
M 



r.(i^' 



H + x; 4) 



(35) 



H 



Transitions involving excited states very near their kinematic threshold for 
an allowed decay (e.g. where the light pseudoscalar momentum is less than 
100 MeV) are extremely sensitive to our calculated mass values. This is 
particularly true for the allowed transitions within the IS" multiplets. In 
these cases, even the small mass differences between the the charged and 
neutral states are important. Using the physical masses for the various D 
0, the individual pion transitions are shown in Table ^ 



mesons 



After removing these phase space uncertainties, reasonable variations in 
our model parameters gave variations of about 10% in the overall hadronic 
widths. The listed branching ratios with emission of a tt are flavor blind 
and sum over the final state pion charge (i.e. they have been computed 
with ( = Vs). Each exclusive decay can be deduced by correcting for this 
factor using Table ^ to determine the relative strength of the charged and 
neutral decays. In addition, a small phase space correction should be included 
appropriate to the slight difference in the masses of the various charge states. 
These rates are shown in Table p for the D and Ds mesons and in Table ^ for 
the B and Bg mesons. For the 1 2 P(i.2) Bg states, our model predicts that they 
are below threshold for K transitions to the the l5S'(o,i) B states. However, 
this is very sensitive to the details of the model. So, for completeness, we 
note the partial rates divided by the appropriate phase space factor at px = 
in Table 0. 

A list of the allowed transitions for other low-lying excited states is re- 
ported in Appendix B. 
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H'{n'^'£j,) H{nHj) x 


4 P. Tj{g\Y 


D(l^Po) ^(l^^o) vr 


437 189 


P'(ltPl) D{l^2Si) TT 

ollhSi) TT 


355 (*) 1.7 
2 355 14.5 


P)(liP2) D{I^Sq) TT 
P)(l55l) TT 


2 506 24.6 
2 394 13.7 


P)(l5Pi) P)(l55i) TT 


420 181 


D,(l5Po) D{lhSo) K 


325 236 


P>,(l5Pi) P)(li5i) 

D{l^Si) K 


175 (*) 1.89 
2 175 0.3 


P>,(liP2) D(l55o) 

D{l'2Si) K 


2 442 8.9 
2 264 1.4 
2 248 0.4 


Ds{l-2Pi) D{l-2Si) K 


302 224 



Table 8: The heavy-light IP state hadronic transition rates for D and Dg 
mesons. H' ^ H + x. Decays denoted with an (*) are allowed only because 
of the order l/m^ mixing of states, p^ and rx/{g\y are in MeV. 



H'{n'^'lj,) H{nHj) x 


4 p. r,/(^8)2 


BiliP,) P(li5„) 7, 


388 186 


P(liPi) B{l^Si) TT 

BillSi) TT 


338 (*) 0.5 
2 338 13.1 


P(15P2) P(l5S'o) TT 
B{l^Si) TT 


2 396 10.6 
2 352 9.5 


P(l5Pi) P(l5Si) TT 


381 180 


P,(l5Po) P(l5S'o) K 


170 159 


p,(i^Pi) p(i^5i) /s: 


153 143 



Table 9: The IP state hadronic transition rates for B and Bg systems. H' — > 
H + X. Decays denoted with an (*) are allowed only because of the order 
l/rrih mixing of states. Values for and Tx/igAY are in MeV. 
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H'{n'^'i,r) H{nHj) x 


4 r,./(^7i)2 X (100M)(2'^+i) 


fi,(liPi) K 


(*)4.92xl0-i 


5,(liPi) B{l-2Si) K 


2 2.38 X 10-2 


5,(liP2) 5(l55o) K 


2 9.48 X 10-3 


P,(liP2) P(1^5i) 


2 1.43 X 10-2 



Table 10: Decay rates for laP^^ 2) Ps mesons with phase space dependence 
divided out. These states are very near the kinematical threshold in our 
model. Values for and Vx/^g^)"^ are in MeV. 

3.3 Comparison to lattice results 

As one more consistency check of our model we compare our prediction for 
the transition! P — P* + vr with model independent results coming from 
lattice simulations. We define 



A 



BB*7 



- E 

3 ^ 

M=l,2,3 



{B*\A,,{v)\B) da 



(36) 



with A^ir) = g(t, r)7^7^g(t, r) at t = 0. The same matrix element can be 
expressed in terms of our radial wavefunction and, in the limit — > 0, 



A 



BB*7 



-9\ 



(37) 



g\ is the effective coupling of the transition as defined in Eq. (p8DP| In the 
chiral quark model g\ is an effective parameter and it is has to be given 
as input. On the other side, in the context of lattice computations, the 
matrix element in Eq. (^) follows directly from first principles (the QCD 
Lagrangian) and can be computed explicitly. By fitting the lattice results of 
Ref. ||2^ with our prediction for A^^*'^{r) as function of g\ we are able to 



''Even if this transition is kinematically forbidden it is physically relevant to the B —> 
t: + £ + D exclusive decay under the assumption of vector meson dominance. 

^Note from Eq. ( |37| ) that in the non relativistic limit the coupling coincides with 
the coupling constant g that appears in the Heavy Meson Chiral Lagrangian The 
lattice result for this quantity is g = 0.42 ± 0.09. 
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Figure 5: Comparison between a prediction of our model and the lattice 
QCD result 



determine]^ 

g\ = ZaqT'''""^ = 0-53 ± 0.11 (38) 
where g^^^^^^^^^^ jg naive lattice result extracted from the fit and Za 



0.78 is lattice matching factor discussed in Ref. [0]. The error includes 
the statistical error due to the simulation and the fits (~ 10%), and an 
estimate of the systematic error in the matching coefficient and in the the 



chiral extrapolation. This is a preliminary result, as those of Ref. are, 
because of the small lattice size and the poor chiral extrapolation. In any 
case our result is in agreement with lattice determination of Ref. , g\ = 
0.61±0.13, and with experimental results from nucleon and hyperon (3 decays, 

= 0.58 ± 0.02, p6[. A more precise lattice determination is possible and 
may be carried out in the near future. 

Apart for the overall normalization given by g% the radial dependence 
of the function AB^Birir) is predicted independently by our model and by 
the lattice computation. In Figure ^ we present a comparison between our 



5 



We set to zero terms of the order ©(ptt) for consistency with the lattice computation. 
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analytical result, with adjusted normalization, and the lattice data. We 
believe this comparison provides a satisfactory consistency check of the two 
methods^. 



3.4 Comparison to experiment 

The total width of the D*^ meson has recently been measured by the CLEO 
Collaboration |2^. They obtain V = 96 ib4stat ± 22syst(keV). Combining this 



measurement with the well-known branching ratios for the various pionic 
transitions, gives a measurement of chiral coupling constant. We obtain 

g\ = 0.82 ± 0.09 (39) 

within our model. 

Within the chiral quark model, the coupling determined from any tran- 
sition should agree, (i.e the coupling is independent of particular initial or 
final heavy-light states.) Table |ll] lists the various determinations of g\ from 
existing data. Within the existing large uncertainties the various determina- 
tions are consistent. 



4 Conclusions 

We have computed the spectrum and hadronic decay width of the excited 
D, Ds, B and Bg mesons using a relativistic quark model for the masses 
and wavefunction of the heavy-light mesons. This work is based on that 
of refs. [|, H, ^ but departs from these previous works because we choose 
a simpler form for the potential and determined its parameters exclusively 
from fitting the experimental heavy-light spectrum. Moreover we computed 
all corrections within the model to order l/ruh, including mixing between 
nearby states with the same J^. 

Our spectrum results agree very well with the existing data in the D and 
Dg systems. The agreement in the B and Bg systems is also fairly good but 
the experimental situation for the P states is not yet completely resolved. 

For example, our model predicts a spin-orbit inversion for the excited P- 



wave states in these systems. This agrees with the recent CLEO fT^ results 



^ It also suggests a possible use of chiral quark model wavefunctions to isolate excited 
states contributions in lattice correlation functions. 
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state 


width (exp) [Ref.] 


V/{g\f (model) 


9\ 


D+{l-2Si) 


96 ± 4 ± 22 (keV) 


m 




143 (keV) [a] 


0.82 ±0.09 


DO(liPi) 
D+(liPi) 


18.91^1 
28 ±8 






16 
16 


-'-•'-'^-0.11 
1 00+0.18 
-'-•'-'^-0.27 


D+(l5Pi) 


2mt\f ± 26 ± 36 


14 




181 


1.27 ±0.22 




23 ±5 






38 


0.77 ±0.08 




25 ±8 






38 




Ds0Pi) 


< 2.3 






2.0 


< 1.07 


Ds{l"^P2) 


15 ±5 






10.9 


1 1 7+0.18 
' -0.11 


P(l^Pi) 
P(l^Pi) 

P(ltP2) 


73 ±44 

10+15 +29 
^°-l2, -23 

41 ±43 


HN 
HIN 
II3N 


180 

14.0 

20.0 


0.64+°iI 

1.431?:^^ 



Table 11: The 15 and IP heavy-light hadronic transition rates compared to 
experiment. All widths are in MeV unless otherwise indicated. Experimental 
values from PDG [|ll|] unless otherwise indicated. 

[a] Theoretical value corrected for phase space observed mass(see Table |^ 
and the 1.6% branching ratio to + 7 |jTl[ . 

[b] Experimental results depend strongly on model dependent assumptions. 
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for the D mesons; but is in disagreement with prehminary L3 results \V1\ for 
the B mesons. However, in our model we find that the splitting between the 
ji = 1/2 P-wave states (J = 0, 1) more than twice as large as the splitting 
of the ji = 3/2 P-wave states (J = 1,2). which is inconsistent with the 
assumption of equal splitting used by L3 in their analysis. 

The information on the spectrum and the wavefunction was used to com- 
pute a complete list of allowed hadronic decays for these excited mesons into 
lower energy states with emission of a light pseudoscalar meson (77, vr or K) 
and their relative branching ratios. 

We also compared our model prediction for {B*\A*^{r)\B) with lattice 
results. We find good agreement between the shapes of the decay amplitudes. 
We used this comparison to extract a preliminary lattice determination of 
g\, the effective coupling of the chiral quark model. We find the value g\ = 
0.53 ± 0.11. Using this value for g\ to translate the total hadronic decay 
widths in physical units, we would conclude that the widths of the D meson 
P waves are consistently below experimental findings. Of course, the lattice 
results have not yet been extrapolated to the continuum. 

Comparison of our results with a recent experiment measurement of the 
D*^ width yields g\ = 0.82 ±0.09. Using this value, we find much better 
overall agreement with experimental findings. 

Finally, there is also a whole set of hadronic decays with emission of light 
vector mesons which we did not include in our study. For the higher excited 
states these transitions give an important contribution to the total physical 
widths. We intend to study these transitions in a future paper. 
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Appendix A 



In order to be able to factorize the expression for the transition amphtude, 
Eq. (|3TD into radial and angular parts, it is convenient to adopt the following 
basis for the F matrices 



poo ^ ( A = ( ° ""M (40) 

.10 _ / \ yll^f ^ ^ 

CJ„ r ^ V a„ 



where (Tq is the 2x2 identity and 0", are the usual Pauli matrices. 

In this basis we obtain simple expressions for the cf''^{X) coefficients^ 



cfj\X) = cotr(7°XC)(/,m',f|a^n,„Jj,m,£) (43) 



cff{X) = cotr(7°XC)(/,m',f|a^n,„Jj,m,2j-£) (44) 



cff{X) = cotr(7°XC)(/,m',2/-£V''n,„Jj,m,£) (45) 



cf:\X) = cotr{^'XTf){f,m',2j'-i'\a^Yk,^Jj,m,2j-i) (46) 



where 



Co = ^A-^Y^ J, J, ^2 ^^^^ 

For the particular case X = jrf' the angular dependence from the p vector 
disappears and we obtain the following explicit expression 

oT^in') = +C2\p\{f,i'\\T^'^\\j,i) (48) 

cZ'\^l') = +cm{f,i'\\Yk\\j,2j-£) (49) 

cf/i^l') = -c,po{f,2f-i'\\Y,\\j,e) (50) 

cffi^l') = -C2bK/,2/-£'||T(f||j,2j-£) (51) 



^In this Appendix we follow the notation of Elbaz fesl where 



[ai...a„] = V(2ai + l)-(2a„ + l) (42) 
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with 



2 



C2 = An{-^n-y'-'^^■^'-''[Jk] \ < U i ^ ) (53) 



2 J' J j \ 
Explicit expressions for the Wigner-Ekkart reduced tensors are 

{i'nY,\\ij) = \l^^^{-y^''-'^^ o){i ^- } (54) 



and 




(55) 



Our expression for the spin structure of the decay amphtudes, reported in 
this Appendix, disagree with Ref. in the overall phase factor. This factor 
does not affect their results but is relevant in case of mixing. 



Appendix B 

In this appendix we list the hadronic transitions H' ^ H + x for the low- 
lying excited states not discussed in Section 3.2. The following Tables list 
the transitions in decending order according to: 

• The flavor of the decaying heavy meson if'; 

• Its mass; 

• The decay channel. 

Transitions with a branching ratio less than 1% are not reported. Some 
decays are marked with an asterisk. These are decays that apparently do 
not conserve the heavy quark spin (|j' — j\ < ix < \j' + j\) but, we remind 
the reader, our initial and final states are physical states and therefore such 
decays become allowed because of mixing effects. 
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2 


685 


2.7 


TT 


2 


818 


62.5 


TT 


1 


540 


1.5 


TT 


1 


506 


5.9 


TT 


3 


470 


4.0 


TT 


1 


444 


3.9 


TT 


2 


255 


5.1 


TT 





138 


3.8 


K 


2 


620 


9.9 



32 



H' = D{2iP2) 



3.035 GeV 



H' = DiliPi) 



H{nHj) 



D{lhSi) 
D{l^Si) 

D{liP2) 
D{liP2) 
D{l--Pi) 

D{liD2) 



X 


n 


Px 




V 


2 


828 


3.3 


TT 


2 


936 


53.4 


V 


2 


721 


2.3 


TT 


2 


848 


47.1 


TT 


3 


541 


4.3 


TT 


1 


505 


4.2 


TT 


3 


505 


2.2 


TT 


1 


480 


6.1 


TT 


2 


393 


10.8 


TT 


2 


296 


5.8 


TT 





142 


3.9 


K 


2 


779 


15.3 


K 


2 


658 


8.7 



H' = D(2^Pi) 



3.045 GeV 



H{nHj) 



D{l2Si) 
d{i^Si) 
D{liPo) 
DiliPi) 
D{liP2) 
D{l--Pi) 
D{2^Si) 



X 


e 

^x 


Px 


^xl{g\f 


V 





730 


10.9 


TT 





855 


85.4 


TT 


1 


582 


61.0 


TT 


1 


549 


6.6 


TT 


1 


513 


19.9 


TT 


1 


488 


61.5 


TT 





306 


50.6 


K 





667 


59.3 



H' = D{UFz) 



H{nHj) X 4 



3.074 GeV 

7mf 



Px 



D{l^Si) 
d{i^Si) 
D{l-2Po) 
D{l"^Pi) 
D{liP2) 
Dil^Pi) 
D{liD2) 
D(liL>i) 
Ds{l--Si) K 



4 
4 
3 
3 
3 
3 
2 
2 
4 



754 
875 
606 
573 
538 
513 
251 
230 
693 



0.5 
13.1 
11.4 
7.3 
3.8 
5.0 
2.9 
3.3 
1.5 



H{nH.j) X l, 



rn = 3.091 GcV 



Px 



D{liSi) 
DiliPi) 
D{liP2) 
D{liPi) 
D{liD2) 
D{liD3) 
D{liD2) 



rj 

TT 
TT 
TT 
TT 
TT 
TT 
TT 
TT 

K 
K 



872 
974 
887 
587 
552 
527 
268 
244 
208 
824 
707 



0.6 
9.8 
7.8 
4.4 
9.9 
9.6 
0.6 
2.9 
2.4 
2.0 
0.9 



H' = D{l^F2) 



m ■ 



3.101 GeV 



H{nHj) 



D{1-^Sq) 

D{liPi) 
D{liPi) 
D{liP^) 
D(liP2) 
D{liP2) 
D{liPi) 
D{liD2) 
D{liD2) 



X 


4 


Px 


^xl{g\Y 


TT 


2 


981 


3.7 


TT 


2 


895 


1.8 


V 


1 


362 


3.7 


TT 


1 


595 


22.2 


TT 


3 


595 


1.4 


TT 


1 


560 


2.2 


TT 


3 


560 


2.1 


TT 


3 


536 


2.6 


TT 





279 


56.8 


TT 


2 


279 


3.4 


K 


2 


832 


3.1 


K 


2 


716 


1.2 


K 


1 


250 


6.0 



H' = D{l-iF3) 



m ■ 



3.123 GeV 



H{nHj) 



D{l2Si) 
D{liPo) 
D{liPi) 
D{liP2) 
D{liP2) 
D{liP2) 
D{l--Pi) 
D{liD2) 

D{liD:i) 
Ds{liP2) 



X 


4 


Px 


Tx/{9\Y 


TT 


2 


910 


4.9 


TT 


3 


645 


4.3 


TT 


3 


613 


3.2 


V 


1 


332 


3.4 


TT 


1 


579 


23.7 


TT 


3 


579 


2.6 


TT 


3 


554 


1.2 


TT 


2 


300 


1.2 


TT 





277 


56.7 


TT 


2 


277 


3.6 


K 


2 


735 


3.3 


K 


1 


203 


3.7 



33 



3.141 GeV 



F'==D,(liDi) 



2.913 GcV 



H{nHj) 


X 




Px 


r./(5l)^ 


D{l^Si) 


V 


1 


811 


1.7 


D{l^Si) 


TT 


1 


923 


44.5 


D{l-2Po) 


TT 





660 


2.9 


D{liP2) 


TT 


2 


594 


15.1 


D{2^Sx) 


TT 


1 


396 


16.3 


D{2^Po) 


TT 





128 


6.5 




K 


1 


751 


8.7 




K 





384 


1.4 



H' = D{^^Si) TO = 3.226 GeV 



H{nHj) 



^(liS'i) 
D{\ip^) 
D{liP2) 

D{2^Si) 
D{2iPi) 



X 




Px 


^x/ig^AY 


1] 


1 


976 


3.2 


TT 


1 


1066 


50.2 


TT 


1 


984 


53.5 


TT 


2 


697 


23.1 


TT 


2 


663 


17.0 


TT 


1 


560 


41.3 


TT 


1 


473 


33.5 


TT 





113 


5.6 


ii" 


1 


930 


17.6 




1 


822 


13.4 



H' = D,{2iSo) 


TO 


= 2.700 GeV 


H{nHj) X 


^x 


Px T^/{9\Y 


D{l^Si) K 


1 
1 


424 3.12 
187 0.04 


H' = D,{2hSx) 


TO 


= 2.806 GeV 


H{nHj) X 




Px T,J{g\Y 


D{l^So) K 
^(liS'i) K 


1 
1 

1 
1 


661 21.1 

539 12.2 

540 6.2 
371 1.5 


H' = D,{liD2) 


TO 


= 2.900 GeV 


H{nHj) X 


^x 


Px T,,/{g\Y 


D{l^Si) K 


3 
3 


629 10.6 
486 1.4 



H{nHj) 


X 




Px 






K 


1 


752 


26.1 


D{l^Si) 


K 


1 


641 


10.7 


D{l^Pi) 


K 





47 


39.7 




V 


1 


644 


15.2 




V 


1 


501 


4.5 



D{l2So 
D{l^Px 
D{l^Pi 



H' = D 



DiliPo 
DiliPi 
D{liP2 
D{liPi 



H' = D, 




TO 


= 2.925 GeV 


H{nHj) 


X 


ix 


Px r,/{g^^r 




K 
K 
) V 
) V 


3 
3 
3 
3 


762 11.4 
652 7.3 
656 3.1 
514 1.1 


H' = Ds 




m 


= 2.953 GcV 


Hin'Cj) 


a: 






D{liSi) 


K 
) l 


1 
1 


677 35.0 
543 16.4 



22P0) m = 3.067 GeV 



X 




Px 


Tx/ig'AY 


K 





875 


74.1 


K 


1 


377 


32.3 


K 


1 


270 


17.8 







780 


49.1 



2iPi) TO = 3.114 GeV 



X 


^X 


Px 


Tx/ig^y 


K 





813 


(*) 0.7 


K 


2 


813 


36.6 


K 


1 


482 


5.5 


K 


1 


436 


1.6 


K 


3 


383 


0.9 


K 


1 


343 


4.5 




2 


700 


11.1 


V 


1 


275 


1.6 



34 



H' = D,{2iP2) m = 3.157 GeV 



H{nHj) 


X 




Px 


rx/(3l)^ 




K 


2 


943 


34.9 




K 


2 


848 


28.9 


D(liPi) 


K 


3 


485 


1.8 




K 


1 


400 


11.0 


£1(2^ So) 


K 


2 


255 


1.3 




V 


2 


854 


14.8 




V 


2 


738 


9.5 


H' = DA 


25 Pi) 


m 


= 3.165 GeV 


H{nHj) 


X 




Px 


^xl{g\Y 


D{l^Si) 


K 





854 


72.1 


D{l-2Po) 


K 


1 


537 


42.6 


D{liPi) 


K 


1 


494 


12.6 


D{llP2) 


K 


1 


446 


41.3 


£)(l5Pi) 


K 


1 


410 


28.7 




V 





745 


45.1 


DsiliPo 


) V 


1 


357 


8.1 


DsiliPi 


V 


1 


280 


2.7 


D,iliP2 


V 


1 


184 


3.4 


H' = DA 


1-2 Fs) 


m 


= 3.203 GeV 




X 




Px 


^xl{9\Y 


D{l^Si) 


K 


4 


884 


9.2 


D{l-2Po) 


K 


3 


575 


5.5 


D{llPi) 


K 


3 


534 


3.2 


D{liP2) 


K 


3 


489 


1.5 


D{liPi) 


K 


3 


456 


1.6 




V 


4 


778 


2.2 


Dsil^Po 


V 


3 


409 


0.5 


H' = Ds( 


15 P4) 


m 


= 3.220 GeV 


Hin^e.j) 


X 




Px 


^xl{9\Y 




K 


4 


990 


7.2 


D{l^Si) 


K 


4 


897 


5.6 


D(liPi) 


K 


3 


552 


2.0 


D(liP2) 


K 


3 


508 


4.0 




K 


3 


475 


3.3 




V 


4 


904 


2.4 




V 


4 


793 


1.4 



H' = Ds{\iF2) 




= 3.224 GcV 


11 1 


X 




Px 


^ x/\9a) 




K 


2 


993 


5.4 




K 


2 


900 


3.1 


P'(liPi) 


K 


1 


556 


39.1 


P>(liPi) 


K 


3 


556 


0.9 


i:'(iiP2) 


K 


1 


512 


4.0 




K 


3 


512 


1.1 


DiliPl) 


K 


3 


480 


0.9 




V 


2 


907 


3.1 




V 


2 


796 


1.5 


D^l^Pi) 
D4liP2) 


V 


1 


372 


14.0 


V 


1 


302 


0.9 



H' = D 



D{l2Si 
D{liPo 
D(liPi 
D{liP2 

D{llP2 

D,iliP2 



H' = D 



D{l2Si 
D{liPo 
D{liP2 



liPg) 



m 



3.247 GeV 



X 


4 


Px 


^x/ig'Xr 


K 


2 


918 


8.2 


K 


3 


619 


2.2 


K 


3 


580 


2.0 


K 


1 


536 


42.2 


K 


3 


536 


1.5 


V 


2 


815 


3.9 


V 


1 


339 


12.6 



3^5o) 



3.259 GeV 



X 


ix 


Px 




K 


1 


927 


19.0 


K 





630 


2.1 


K 


2 


549 


5.1 


K 


1 


254 


1.0 


V 


1 


825 


6.8 


n 





478 


U.9 



35 



3.345 GeV 



H' = B{liDi) 



n \ii tj ) 


X 




Px 


r Kn^V 
J- a:/ 15a j 




K 


1 


1080 


28.0 


D{liSi) 


K 


1 


993 


26.6 


-0(1 5 Pi) 


K 


2 


675 


12.2 


K 


2 


635 


8.1 


D{llPi) 


K 





607 


1.8 




K 


1 


507 


18.8 


D{2^Si) 


K 


1 


385 


8.4 




V 


1 


1001 


13.6 




V 


1 


896 


10.6 


Ds{liPi) 


V 


2 


523 


1.7 


Ds{2iSo) 


V 


1 


309 


1.7 



H' 



Bil^So 
B{liSi 
5(1 i Pi 
5(1 iP2 
5(1 5 Pi 



H' = 5( 


2^ So) 


m 


= 5.8 


86 GeV 


H{nHj) 


X 


4 


Px 


Tx/ig^r 


5(l^S'i) 
5(l5Po) 


TT 
TT 


1 




519 
114 


21.9 
4.9 


H' = 5( 


2^5i) 


m 


= 5.920 GeV 


H{nHj) 


X 


4 


Px 


r./(g^)^ 


B{l^So) 
bIi'^Si) 
5(1 i Pi) 

5(l5Pi) 


TT 
TT 
TT 
TT 


1 
1 

2 



591 
550 
167 
109 


19.2 
22.9 
0.5 
4.6 


H' = 5( 




TO 


= 5.985 GeV 


H{n^(.j) 


X 


4 


Px 


rx/ig^kr 


B{l^Si) 
5(1 i Pi) 
5(l5Po) 
5(liP2) 

5(l5Pi) 


TT 
TT 
TT 
TT 
IT 


3 
2 
2 
2 
2 


611 
244 
237 
228 
195 


17.4 
1.8 
1.6 
0.7 
0.5 





m 


= 5.993 GcV 


X 


f 

^x 


Px 


r./(5i)^ 


TT 


3 


658 


11.1 


TT 


3 


618 


10.6 


TT 


2 


252 


0.6 


TT 


2 


237 


2.2 


TT 


2 


204 


1.3 



H{nHj) X 4 



6.025 GcV 

^^l{g\f 



P. 



B{l2So) 
B{liSo) 
B{llSi) 
B{liPi) 
5(liPi) 
Bsil^So) 
S,(1^5i) 



V 

TT 
TT 
TT 
TT 

K 
K 



475 
687 
647 
286 
286 
403 
330 



2.8 
18.2 
7.5 
81.4 
1.4 
7.2 
2.0 



H' = B{liD2) 



6.037 GeV 



HinHj) 



5(l55i) 
5(1551) 
B{liPi) 
5(liP2) 
5(liP2) 
5,(1^5i) 



X 


f 

^x 


Px 


^x/ig^r 


V 


1 


431 


3.3 


TT 


1 


658 


23.9 


TT 


2 


299 


1.3 


TT 





284 


81.1 


TT 


2 


284 


2.0 


K 


1 


350 


7.1 



5' = 5(22Po) TO = 6.163 GeV 



5(l55o) 
5(l5So) 
B(llPi) 
5(l^Pi) 
5(25So) 



X 




Px 


^x/{g\Y 







643 


9.9 


TT 





810 


97.4 


TT 


1 


425 


21.1 


TT 


1 


383 


52.9 


TT 





233 


39.6 


K 





576 


51.0 



H' = 5(23Pi) 



TO 



6.175 GeV 



B{l2Si) 

5(1551) 

S(liPi) 
5(l5Po) 
5(liP2) 
5(l^Pi) 
5(2551) 
5(li5»i) 

S.(l35l) 



X 


4 


Px 


^x/{g\Y 


V 


2 


606 


1.5 


TT 


2 


782 


59.6 


TT 


1 


437 


1.6 


TT 


1 


431 


2.4 


TT 


3 


423 


2.6 


TT 


1 


395 


3.4 


TT 


2 


209 


2.3 


TT 





55 


1.4 


K 


2 


534 


4.6 



36 



H' = B{2iP2 



m 



6.188 GeV 



H' = B{liF4 



B{liSi) 
B{liSi) 
BiliPi) 
B{liP2) 
B{liPi) 
B{2iSo) 
B{2^Si) 
B{liD2) 
Bsil^So) 
B^il^Si) 



X 


n 


Px 




V 


2 


671 


1.2 




2 


832 


36.0 


V 


2 


621 


1.1 


TT 


2 


793 


39.7 


TT 


3 


449 


2.2 


TT 


3 


436 


1.2 


TT 


1 


408 


5.7 


TT 


2 


261 


2.4 


TT 


2 


224 


1.9 


TT 





57 


1.5 


K 


2 


605 


4.2 


K 


2 


550 


3.4 



ir = B(2^ri_) ;// = (i.l!)4 GcV 



B{l2Si) 
B{l^Si) 
B(liPi) 

S(l3Po) 
B(15P2) 
B(l3Pi) 

B{2^Si) 
Bs{l^Sx) 



V 

TT 

TT 
TT 
TT 
TT 
TT 

K 



629 
799 
455 
449 
442 
414 
231 
557 



9.2 
93.2 

6.1 
30.6 
22.6 
39.5 
38.7 
46.3 



H' = B{l2F3) 



6.220 GeV 



H{nHj) 



B{l2Si) 
B{liSi) 
Bil-^Pi) 
B{l^Po) 
B{liP2) 
B(l3Fi) 
B{liD2) 
B{liDi) 
Bs{l^S^) 



X 


f 


Px 




V 


4 


659 


0.3 


TT 


4 


822 


11.7 


TT 


3 


481 


3.5 


TT 


3 


475 


3.8 


TT 


3 


468 


1.9 


TT 


3 


440 


2.1 


TT 


2 


186 


0.8 


TT 


2 


135 


0.3 


K 


4 


588 


0.5 



6.226 GcV 



P. 



B{l^So) 
B{l^Si) 
B(liPi) 

5(1^^2) 

B{liPi) 
Bil^Ds) 
B,{liSo) 



TT 
TT 
TT 
TT 

TT 
TT 

K 
K 



4 
4 
3 
3 
3 
2 
4 
4 



865 
827 
486 
473 
445 
183 
647 
594 



7.0 
6.7 
1.5 
4.7 
4.4 
0.8 
0.5 
0.3 



H' = B{liF2) 



m ■ 



6.264 GeV 



H{nHj 



B(l5^o) 
B{lhSi) 
B{liPr) 
S(liP2) 
B{riP2) 
B(l^Pi) 
B{liD2) 
B{\iD2) 
Bs{l--So) 
B,(1^5i) 



X 


4 


Px 




TT 


2 


898 


3.7 


TT 


2 


860 


2.0 


TT 


1 


522 


18.9 


TT 


1 


509 


2.0 


TT 


3 


509 


1.3 


TT 


3 


482 


1.3 


TT 





236 


51.2 


TT 


2 


236 


1.6 


K 


2 


688 


1.9 


K 


2 


637 


0.9 



H' = B{liF^) 



6.271 GeV 



HjnHj) 



B{l^Si) 
S(l55i) 
B(l^Pi) 
S(l^Po) 
B{liP2) 
B{liP2) 
B{1^D3) 
Bil^Ds) 
Bsjl^Si) 



X 


ix 


Px 




V 


2 


714 


1.0 


TT 


2 


866 


5.2 


TT 


3 


529 


1.2 


TT 


3 


523 


1.0 


TT 


1 


516 


20.4 


TT 


3 


516 


1.5 


TT 





235 


51.0 


TT 


2 


235 


1.7 


K 


2 


644 


2.3 



37 



H' = B{3iSo 



B{l2Si) 
B{\^Si) 
B{l-2Po) 
S(liP2) 
5(2^51) 
B{2^Po) 
Bsil^Si) 



So) 


m 


= 6.320 GeV 


X 




Px 




V 


1 


766 


1.1 


TT 


1 


908 


49.5 


TT 





569 


1.8 


TT 


2 


562 


11.0 


TT 


1 


363 


9.5 


TT 





72 


3.5 


K 


1 


697 


5.1 



H' = B,{liD^) 



6.103 GcV 



H' = B{3iSi) 



H{nHj) 



6.347 GeV 

Px T,/{g\Y 



B{l2So) 
B{lhSi) 
B{liPi) 

B{llP2) 
S(l3Pi) 

B{2iSo) 
B{2^Si) 
B{2^Pi) 
BsiliSo) K 
B.(l55i) K 



970 

932 
600 
587 
560 
423 
390 
64 
775 
726 



31.3 
43.0 
9.4 
7.9 
1.7 
9.3 
10.8 
3.0 
5.0 
5.2 



ff' = B,,(2^5o) 


m 


= 5.985 GeV 


H{n-'f'j) X 






B{l^Si) K 


1 


41() 3.15 


H' = B,{2iSi) 


m 


= 6.019 GeV 


H{n^ej) X 


'^x 


Px T,/{g\Y 


B{liSo) K 
B{\^Si) K 

S.(1^5o) V 


1 
1 
1 


517 5.6 

462 4.9 
325 0.4 


H' = B.s{liD2) 


m 


= 6.095 GeV 


H{n3£j) X 




Px T,/{g\Y 


B{l^Si) K 
BsiliSi) rj 


3 
3 


555 6.62 
370 0.30 



H{n^lj) 


X 




Px Vxl{g\Y 


B{l^So) 
5(155-1) 

Bsil^So 
Bs{l^Si 


K 

K 
V 
V 


3 
3 
3 
3 


614 5.03 
564 4.17 
453 0.46 
383 0.21 


H' = Bs( 


l*i?i) 


m 


= 6.127 GeV 


H{nHj) 


X 


^x 


Px ^xl{g\Y 


5(1 5 5o) 

5(15 ^l) 

5,(15^0 
5, (15^1 


K 

K 
V 
V 


1 
1 
1 
1 


641 26.9 

592 11.5 
486 10.1 
420 3.4 


H' = BA 


lifz) 


m 


= 6.140 GeV 


H{nHj) 


X 




Px ^xl{g\Y 


5(1 2 Si) 
5,(l35i 


K 
V 


1 
1 


607 36.2 
438 11.5 


H' = 5,( 


2iPo) 


m 


= 6.264 GeV 


H{nHj) 


X 


ix 


Px rj(5^)^ 


5(1550) 
5(1 iPi) 
5(l55i) 
5,(1^50 


K 
K 
K 
V 



1 
1 



785 77.7 
262 20.6 
164 3.8 
652 40.5 



n' = B. 



5(l25i 
5(l55i 

5(l^Po 

5(l55i 

5,(l55i 



H' = B, 



5(l55o 
5(l55i 
5(1 5 Pi 

5.(l5 5o 
5,(1551 



25Pi) 



6.278 GeV 



X 


ix 


Px 


rx/ig^Ar 


K 





755 


(*) 0.7 


K 


2 


755 


28.7 


K 


1 


278 


3.3 


K 


1 


201 


1.4 


V 


2 


614 


5.8 



2 2 Pa) TO = 6.292 GeV 



X 


ex 


Px 


rx/{g\Y 


K 


2 


813 


19.2 


K 


2 


769 


19.8 


K 


1 


232 


4.2 


V 


2 


683 


4.8 


V 


2 


630 


4.2 
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H' = B. 



B{liPi 

B{liP2 
B{l^Pi 
Bsil^S, 



H' = B. 



B{l^Si 
5(1 i Pi 

5(1 iP2 
5(1 5 Pi 

5,(l3 5l 



H' = B, 



B{l2So 
B{liSi 
5(1 i Pi 
5(1 iP2 
5(1 5 Pi 

5,(1^^0 



H' = 5; 



H{nH„ 



B{l2So 
B{l^Si 
B(liPi 
5(liP2 
5,(l5 5o 
5,(l5 5i 
S,(liPi 



2^ Pi) 



6.296 GeV 



X 




Px 


r./(5l)^ 


K 





773 


73.8 


K 


1 


319 


7.0 


K 


1 


309 


9.3 


K 


1 


296 


24.5 


K 


1 


241 


6.0 


■n 





634 


37.1 



liPg) m = 6.332 GeV 



X 




Px 


^xl{g\Y 


K 


4 


808 


7.07 


K 


3 


375 


0.51 


K 


3 


366 


0.50 


K 


3 


354 


0.22 


K 


3 


307 


0.14 


V 


4 


674 


1.04 





m 


= 6.337 GeV 


a: 


(.r 


P.r 




K 


4 


857 


4.5 


K 


4 


814 


4.1 


K 


3 


383 


0.2 


K 


3 


363 


0.6 


K 


3 


317 


0.3 


) V 


4 


732 


0.8 


) V 


4 


681 


0.6 



1*P2) 



m ■ 



6.369 GeV 



X 


4 


Px 


^x/ig^AY 


K 


2 


888 


6.0 


K 


2 


845 


3.5 


K 


1 


428 


29.6 


K 


1 


409 


3.0 




2 


766 


2.4 


V 


2 


715 


1.3 


V 


1 


129 


0.9 



H' = B, 


itPa) 


m 


= 6.376 GeV 


H{nHj) 


X 


^x 


Px r^/igW 


B{l^Si) 


K 


2 


853 9.0 


5(liP2) 


K 


1 


420 31.6 




V 


2 


723 3.3 



ff' = 5,(3^^o) 



6.421 GcV 



P. 



5(125-1) K 

Bil^Po) K 

5(liP2) K 

Bs{l--S^) Tj 



1 896 

489 

2 479 

1 770 



16.5 
0.3 
1.7 
3.9 



F' 5,(3^S'i) m = 6.449 GeV 



H{nHj) 



B{l^So) 
P(l'5i) 
B{liPi) 
5(liP2) 

Ps(l^^o) 

Bs{l^Si) 



X 




Px 


^^/{gW 


K 


1 


963 


13.2 


K 


1 


922 


16.0 


K 


2 


531 


2.2 


K 


2 


514 


1.7 


V 


1 


848 


3.9 


V 


1 


800 


4.1 
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